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ABSTRACT 

Three energy mechanisms invoking large-scale magnetic fields are incorporated in 
a model to interpret jet production in black hole (BH) systems, i.e., the Blandford- 
Znajek (BZ) , the magnetic coupling (MC) and Blandford-Panye (BP) processes. These 
energy mechanisms can coexist in BH accretion disc based on the magnetic field con- 
figurations constrained by the screw instability, provided that the BH spin and the 
power-law index indicating the variation of the magnetic field at an accretion disc are 
greater than some critical values. In this model the jets are driven by the BZ process 
in the Poynting flux regime and by the BP process in the hydromagnetic regime, being 
consistent with the spine/sheath jet structure observed in BH sources of stellar and 
supermassive size. 

Key v^rords: accretion, accretion discs — black hole physics — magnetic field — jet 
production 



1 INTRODUCTION 

, As is well known, jets exist in many astronomical cases, such 
1 as active galactic nuclei, quasars and young stellar objects. 
' Different theoretical models have been proposed for acceler- 
ation and collimation of jets, which can be divided into two 
main regimes. Energy and angular momentum are carried by 
both the electromagnetic field and the kinetic flux of mat- 
ter in the hydromagnetic regime, and those are carried pre- 
dominantly by the electromagnetic field in the Poynting flux 
regime (Ustyugova et al. 2000; Lovelace et al. 2002). Bland- 
ford & Znajek (1977) proposed firstly that jets from AGNs 
can be powered by a rotating black hole (BH) with a large 
scale magnetic field threading its horizon. Later, Blandford 
& Payne (1982, hereafter BP82) suggested that an outfiow 
of matter can be driven centrifugally by large-scale mag- 
netic fields anchored at the disc surface. These two mech- 
anisms are usually referred to the BZ and BP processes, 
respectively. As argued in BP82, an outfiow of matter can 
be driven centrifugally from the disc, provided that the an- 
gle of the poloidal magnetic field with the normal of disc 
surface is greater than a critical value, i.e., a^^ > 30''. The 
BZ and BP processes belong to the Poynting fiux and hy- 
dromagnetic regimes, respectively. 

Recently, much attention has been paid to the mag- 
netic coupling (MC) of a rotating BH with its surround- 
ing accretion disc, and this mechanism is referred to as the 
MC process, which can be regarded as a variant of the BZ 



process (Blandford 1999; Li 2000, 2002a; Wang, Xiao & Lei 
2002, Wang et al. 2003, hereafter W02, W03; Uzdensky 2004, 
2005). In the MC process energy and angular momentum are 
transferred from a rotating BH to its surrounding disc. Al- 
though the MC process cannot power jet /outflow directly, 
it plays an important role in depressing disc accretion due 
to transfer of angular momentum from a rotating BH to the 
inner disc. 

In this paper we incorporate the BZ, MC and BP pro- 
cesses into a model to explain the jets from AGNs and the 
BH binaries. Henceforth this model is referred to as the BZ- 
MC-BP model. It turns out that the three mechanisms can 
coexist, provided that the BH spin and the power- law index 
indicating the variation of the magnetic field with the disc 
radius are greater than some critical values. Since the jets 
could be driven in two regimes (the BZ process in the Poynt- 
ing fiux regime and the BP process in the hydromagnetic 
regime), this model could be applicable to a spine/sheath jet 
structure observed in BH sources of stellar and supermas- 
sive size (Meier 2003). This paper is organized as follows. In 
Sect. 2 the magnetic field configuration of the BZ-MC-BP 
model is outlined. The magnetic field configuration is given 
based on the constraint of the screw instability, which con- 
sists of three parts: region I, II and HI corresponding to the 
BZ, BP and MC processes, respectively. In Sect. 3 the con- 
dition for a centrifugally driven outfiow of matter from the 
disc is discussed in a parameter space. The matter outfiow 
rate is determined based on the radial variation of the accre- 
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tion rate due to the MC effect, and the BP power is derived 
based on the work done by the magnetic torque given in 
BP82. It is shown that the outer boundary radius of region 
II for the BP process is intimately related to the power-law 
index indicating the variation of the accretion rate in region 
III. In Sect. 4 we compare the relative importance of the 
BZ and BP processes in jet production, and estimate the jet 
powers as the sum of the BP and the BZ powers. Finally, 
in Sect. 5, we discuss the potential application of the BZ- 
MC-BP model to astrophysics. It turns out that the main 
features of this model are consistent with the general rel- 
ativistic MHD simulations of accretion and outflow in BH 
systems. 

Throughout this paper the geometric units G = c = 1 
are used. 



2 MAGNETIC FIELD CONFIGURATION OF 
BZ-MC-BP MODEL 

In order to discuss the magnetic extraction of energy and 
angular momentum from BH accretion disc by virtue of the 
BZ, MC and BP mechanisms, we must determine the cor- 
responding magnetic field configurations. The difficulty re- 
lated to this issue lies in two aspects. First, we have not 
enough knowledge about the origin of magnetic field in the 
BH accretion disc. Some authors argued that magnetic field 
could be amplified by virtue of dynamo process at accretion 
disc, and the magnetic field at the BH horizon is brought 
from disc accretion and kept by the magnetic pressure of 
the surrounding disc (MacDonald & Thorne 1982; Ghosh & 
Abramowicz 1997; Balbus & Hawley 1998). Li (2002b) and 
Wang et al. (2007) discussed the origin of the magnetic field 
configuration corresponding to the MC process by assuming 
a toroidal electric current fiowing in the equatorial plane of a 
Kerr BH. Unfortunately, a reliable origin of large-scale mag- 
netic fields in BH accretion discs remains unclear. Second, 
we cannot determine the connection between the magnetic 
field configurations corresponding to the BZ, MC and BP 
mechanisms. 

Lovelace, Romanova & Bisnovatyi-Kogan (1995, here- 
after L95) investigated disc accretion of matter on to a ro- 
tating star with an aligned dipole magnetic field, and they 
argued that when the angular velocities of the star and disc 
differ substantially, the magnetic field linking the star and 
disc rapidly inflates to give regions of open field lines ex- 
tending from the polar caps of the star and from the disc. 
The open fleld line region of the disc leads to the possibility 
of magnetically driven outflows. 

Wang, Lei & Ye (2006, hereafter W06) proposed a 
model to explain the light curves of gamma-ray bursts by 
considering the effects of screw instability of magnetic field. 
It turns out that the screw instability in the BZ process 
(henceforth SIBZ) can coexist with the screw instability in 
the MC process (henceforth SIMC), provided that three pa- 
rameters are greater than some critical values. These pa- 
rameters are (i) the BH spin defined as a* = J/M^, (u) the 
power-law index n indicating the variation of the poloidal 
magnetic field at the disc, oc r^" and (iii) the critical 
height of the astrophysical load He above the disc surface. 
The involved quantities M, J, rd and are the BH mass, 
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Figure 1. Schematic drawing of the magnetic field configuration 
of the BZ-MC-BP model 



the BH angular momentum, the disc radius and the poloidal 
magnetic field at the disc, respectively. 

In this paper we combine the BZ, MC and BP pro- 
cesses into a model by virtue of the screw instability of the 
magnetic field. Based on the constraints due to SIMC and 
SIBZ given in W06 we have the magnetic field configuration 
related to the BZ-MC-BP model as shown in Figure 1. 

In Figure 1 r^^^^^ is the critical radius constrained by 
SIMC, which is determined by 

{2i.t,,,,JLmc)BI/bJ = 1. (1) 

Similarly, r^^^^ is the critical radius constrained by SIBZ, 
which is determined by 

{2i.r,^,JLBz)Bl/Bl^l. (2) 

Equations ([1} and ((2]) are derived based on the Kruskal- 
Shafranov criterion: the screw instability will occur, if the 
magnetic field line turns around itself about once (Kadomt- 
sev 1966; Bateman 1978). In equation ^ Lmc is the critical 
length of the poloidal field line for SIMC, and B^ and Bj are 
the poloidal and toroidal components of the magnetic field 
on the disc, respectively. In equation ([2]) Lbz is the critical 
length of the poloidal field line for SIBZ, and B£ and Bj^ 
are the poloidal and toroidal components of the magnetic 
field on the astrophysical load, respectively. 

As shown in Figure 1, the regions of open field lines at 
the horizon and disc are referred to as regions I and II, which 
correspond to the BZ and BP processes, respectively. The 
region of closed field lines connecting the BH with the disc 
is referred to as regions III, which corresponds to the MC 
process. Region I is confined to the angular region at the 
horizon, < 6 < 0$, while region II is confined to the radial 
region at the disc, r^j^j^. < r < rout- Regions III is confined 
to 9$ < < Ol at the horizon and to Vms < r < r^^j^j^ 
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at the disc. The angle 6$ is the angular boundary between 
the open and closed field lines on the horizon, and 9l is the 
lower boundary angle for the closed field lines. Throughout 
this paper 6l ~ 0.457r is taken in calculations. 

From Figure 1 we find that the magnetic field config- 
uration for the BZ-MC-BP model looks similar to that of 
L95 (see Figure 3 given in L95). In both cases the open field 
lines extend from the central object and from the disc, and 
closed field lines connect the central object with the inner 
disc. However, there are several differences between the two 
cases. 

(1) In L95 the central object is a neutron star, and 
the magnetic field lines are frozen at its surface. The open 
field lines are produced by the toroidal magnetic flux which 
is generated out of the closed fleld lines, arising from the 
difference of angular velocities between the neutron star and 
the disc. While the central object is a rotating BH in the 
BZ-MC-BP model, and the fleld lines can slip on the BH 
horizon. The coexistence of the open and the closed fleld 
lines at the BH horizon has been argued in W06. In the BZ- 
MC-BP model we assume that the open magnetic fleld at 
region II is amplified by the dynamo process, being brought 
inwards by the disc accretion, and the critical radii {rgjj^jf-. 



and Tc 



provide a natural constraint to the open fleld 



lines extending from the disc to inflnity as shown in Figure 
1. 

(2) The directions of the open fleld lines across the sep- 
aratrix are opposite in L95, implying a current sheet exists 
along the separatrix. A magnetic reconnection might start 
from this conflguration. In the BZ-MC-BP model, as shown 
in Figure 1, the open fleld lines in region I and II are in the 
same direction, which can be continuous across the separa- 
trix, and the open magnetic field at region II is balanced by 
the magnetic pressure in the regions I and HI. 

(3) In L95 the magnetic fleld lines in region I penetrate 
the disc vertically, and an outflow cannot be driven centrifu- 
gally in the BP process. Contrary to L95, as shown in Figure 
1, the poloidal magnetic fleld in region II could make an an- 
gle greater than 30*^ with the normal of the disc surface, 
and an outflow of matter driven centrifugally by the open 
magnetic fleld is permitted. 

Following BP82, we assume that the poloidal magnetic 
fleld in region II varies with the disc radius ra as follows, 

- ^-5/4 



(3) 



Considering the balance of the magnetic pressure across the 
boundary of region II and HI, we have 

(Bf = Bmc (r,/r,,„ J-^/^ = Bmc {U^simc)'"^' , (4) 

where Bmc is the magnetic fleld at r^^j^^^, and ^ = rd/vms 
is the disc radius in terms of rms- 

According to W03 the magnetic fleld -Ba/c is related to 
the magnetic field at the horizon, Bh, by 

BH2n {-ojp),^,^ de 



~BMc2n {ujp/VA) ^^^^2 dra. (5) 



where zu, p and A are the Kerr metric coefficients, and they 
read 



zu — (S/p) sin 0, 



E2= (r^ + aY -a^'AsinH, 
p'^ = r^ + cos^ e, 
A = + - 2Mr. 



(6) 



The angular coordinate 6 on the BH horizon is related 
to the dimensionless disc radius ^ by the following mapping 
relation: 



cosy — cos I G(a»;^,n)d5, 



(7) 



where the function G {a,; £,,n) has been given in W03 as 
follows. 



G(a.;5,n) = 



(8) 



In equation (|8]) Xms is deflned in terms of the radius of 
the innermost stable circular orbit, rms ~ Mxms) which is 
regarded as the inner edge of the disc. Incorporating equa- 
tions Q and ((Tj), we have 



Bmc ~ 



Be 



FG(a.;5,n) 



(9) 



where q = a/I — of is a function of the BH spin a,. Thus 
we can determine the poloidal magnetic field in region II by 
combining equation Q with equation (|9]). 



3 OUTFLOW RATE, POWER AND TORQUE 
IN BP PROCESS 

3.1 Condition for outflow driven in 
hydromagnetic regime 

As argued in BP82, an outflow of matter can be launched 
centrifugally, only if the angle of the poloidal component of 
the magnetic fleld with the normal of the disc is greater than 
a critical value, i.e., apL > 30". We can discuss the condition 
for driving the outflow based on the magnetic configurations 
given in Figure 1, and the angle apL can be estimated by 

apL = tan~^ [{^sibz - ^simc) Xms/ he] , (10) 

where he = Hc/M is the dimensionless height of the load 
above the disc surface. 

Incorporating equation (jlOl) with the criterions ([l]) and 
©, we obtain the values of ufl corresponding to the con- 
cerned parameters of SIBZ and SIMC as shown in Table 1, 
in which we have apL > 30° for driving the outflow in the 
hydromagnetic regime. 

It is shown in Table 1 that the angle apL can be de- 
termined by three parameters, a*, n and he, and the corre- 
lation of apL with the three parameters is depicted in the 
a-t — n parameter space for the given he in Figure 2. The 
coexistence of SIBZ and SIMC is indicated by the shaded 
regions, which are bounded by two contours, i.e., 6s = Q and 
he =52.41, 41.55 and 33.73 in Figure 2a, 2b and 2c, respec- 
tively. Each shaded region is further divided by the contour 
of apL = 30°, below and above which we have qfl < 30° 
and apL > 30° in regions A and B, respectively. 

Thus the two regimes of outflows driven can be deter- 
mined by regions A and B, which correspond to the Poynting 
flux regime and hydromagnetic regime, respectively. In this 
paper we confine the discussion to the BP process by taking 
the values of a, and n in region B. 



4 



Table 1. The angle apj^ corresponding to the parameters of SIBZ and SIMC 



a* 


n 


he 




?S/MC 




0.40 


7.6 


52.41 


12.67 


2.20 


42.7° 


0.50 


7.1 


41.55 


10.64 


2.21 


40.7° 


0.60 


6.6 


33.73 


9.31 


2.22 


38.8° 


0.70 


6.1 


27.40 


8.39 


2.23 


37.3° 


0.80 


5.6 


21.55 


7.74 


2.22 


36.7° 




0.4 0.5 0.6 0.7 0.8 0.9 0.4 0.5 0.6 0.7 0.8 0.9 0.4 0.5 0.6 0.7 0.8 0.9 

a«, a« 



(a) (b) (c) 

Figure 2. The shaded regions between the contour of 0s=O for (solid line) and the contours of /ic(dashed lines) are divided by the 
contour of apL=30^ (thick solid line) for he = 52.41, 41.55 and 33.73 in Figures 2a, 2b and 2c, respectively. Regions A and B indicate 
the parameters for a less and greater than 30°, respectively. 



3.2 Rate of matter outflow and BP power and 
torque 

Since angular momentum is transferred magnetically from 
a rotating BH to the inner disc in the MC process, the ac- 
cretion rate is probably depressed. Thus we assume that the 
accretion rate in region II obeys the following relation, 

Md/drd > 0, for r^^^^^ < < rout- (H) 

Based on mass conservation the outflow rate is given by 

Moutfio^^{M,)^^^~{Ma)^,^^, (12) 

where (Md) „ and (Add) are the accretion rates at 

\ / AlC V / out 

'^siMc 8.nd rout, respectively. 

As the magnetic field on the BH horizon is supported 
by the surrounding disc, there are some relations between 
Bh and A'ld- One of them is given by considering the bal- 
ance between the magnetic pressure on the horizon and the 
ram pressure of the innermost parts of an accretion flow 
(Moderski, Sikora & Lasota 1997), i.e., 

Bl/{8n) = Pram ~ pc^ ~ Md/{inrjj), (13) 

As a simple analysis, we take the accretion rate in region 
III as a constant, which is related to the magnetic field at 
BH horizon by the following relation, 

(Md) ^ amBjtrl (14) 



In equation (I14|) the parameter Om is a adjustable parameter 
due to the uncertainty of equation (|13[) . and we take Um ~ 
0.1 in calculations. 

Assuming that the accretion rate in region II varies with 
disc radius in a power-law, Md oc rf , we obtain its expression 
by combining equation (I14p as follows, 

Md = a^B^rlf (rd/rsji^cf /'"' '"s/mc < ^d < rout{l5) 
or 

rhd = Md/BjjM" = am{l + qf {i/£.siMc f 

for isiMC < C < iout 

where BjjM'^ = (Bjmjj x 7.32 x 10^) g- s~^. 

Thus the outflow rate from a ring with width Vd — rd + 
drd can be expressed as 

dmoutfiom = o:,n (1 + q)^ 5* {^/^siMc)^~^ d (^/Cs/A/c) , (17) 

and the total outflow rate driven from region II is 

rhoutfiow = o:,n (1 + g)^ ("jfut - l) , (18) 

where rhoutfiow = Moutfiow/ BjjM^ and ?o„t = ^out/^siMC- 
According to BP82 the specific energy e and angular 
momentum I are constants along each open field line, and 
they read 

e = v^/2 + h + <i>-ujrB^/k, l = rv^-rB^/k. (19) 

The quantities involved in equation (|19p are interpreted 
as follows. The quantity r is the cylindrical radius, and uj is 
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the angular velocity of the foot point of the field line, being 
taken as Keplerian angular velocity of a thin disc and re- 
maining constant along the field line. The quantities h and "I> 
are respectively the enthalpy per unit mass and gravitational 
potential, and the quantities v and are the velocity of the 
streaming gas and its toroidal component, respectively. The 
parameter k is defined as 



(20) 



which is interpreted as the ratio of the mass fiux to the 
magnetic flux, and it also remains constant along each field 
line. The quantities pm and vp are the mass density and 
the poloidal velocity of the streaming gas, respectively. The 
toroidal velocity is related to the toroidal magnetic field 
by 



= vpB^p/Bp + rcu, 



(21) 



where vpB^p/Bp is the toroidal velocity of the streaming 
gas with respect to the rotating field line, and rui is the 
toroidal velocity of the field line itself. Incorporating equa- 
tions (fT9)) — ((2T]) . we have 



B^o = 



kiruo-l/r) _ fcrtj (l - ArJ/r^) 



{vp /ve 



[vp/vAp) 



(22) 



where vap = Bp j -^Airp is the poloidal component of the 
Alfven velocity, and A = ljr\u) is defined as the ratio of the 
specific angular momentum to that at the midplane of the 
disc in BP82. 

Incorporating equations (|19l) — (|22|) . we have 



1-f 



{vplvApf{\~\rllT^) 

1 — (vp jvAp'f' 



(23) 



According to BP82 the term —uorB^/k in equation (|19p 
represents the work done on the streaming gas by the mag- 
netic torque. Combining equation (|22p with the work of the 
magnetic torque and the outfiow rate, we derive the BP 
power and torque as follows. 



dPsp = (-ujrB^/k) dMoutfiow 

2/, 2 2\ 
uj IXr^—r j , 

~ rdMoutfiow, 



dTpp — dPsp/io — 



(vp/vAp) 



;dMo 



utflow • 



(24) 



(25) 



In order to avoid the infinite BP power and torque at 
the Alfven surface as up ^ vap, we have 

A = (rA/raf ■ (26) 

And equations (|24l) and (|25l) can be rewritten as 



dPpp = 



2/2 



(vp/vAp) 



■2dMoutflou 



dTpp — dPpp/uJ — 



(vp/vAp) 



rdM, 



outflow • 



(27) 



(28) 



Since the velocity of the gas increases with the cylinder 
radius in approaching the Alfven velocity, we assume that 
Vp varies with r as follows. 



vp/vAP = (r/r-A)" ; 



(29) 



where a is a parameter to be determined. Substituting equa- 
tion (29) into equation (27) and taking limit for r — > by 
using L'Hospital law, we have 



dPpp ~ lim 



:dMoutfiow ~ dMoutfiow{-iO) 



r^rA 1 — (r/r^) 

Considering the fact that B^p is dominative over Bp 
near the Alfven surface, we infer that vp « v^ and take 
the specific kinetic energy of the streaming gas as LiP'r^j a 
with a = 2 in equation (30) . Substituting equation (26) into 
equation (30), we express the BP power and torque at the 
Alfven surface as follows. 



dPpp = (A/2) a; r^dA^outftou 
dTpp — {X/2) ujr'^dMoutfiow 



(31) 
(32) 



Incorporating equation (17) and integrating equations 
(31) and (32) from ^simc to £,out, we have 



= (Aq„/2) 



S-1 

s 



( S+1/2 
\ ^out 



Tbp = (Aa™/2) (1 + qf Xm4s/MS S + 1/2 

where Pbp = Pbp/Po, Tbp = Tbp/Tq and 



(33) 



(34) 



Po = (B 



p ^2 



AP 



BlruH X 6.59 x 10^® erg ■ 



To = {Bjj) ^ Blm% x 3.26 x lO^^g ■ cm^ • ''^^^ 

The powers and torques in the BZ and MC processes 
have been derived in W02 as follows. 



Pbz = Pbz/Po = 2ai 



fc(i - fc)sin^ ede 



(1 — q) sin^ 



Tbz = Tbz/To = 4a. (1 + q) 



' (l-/c)sin^ 



2-(l-g) sin^e* 



P -P /P / /3(l-/3)sin^ 

Pmc = Pmc/Po = 2a, 



Tmc = Tmc/To = 4a, (1 -f g) 



(1 - g) sin^6i 

(1-/3) sin^ ede 
2- (1-g) sin^e 



(36) 
(37) 
(38) 
(39) 



where k and /3 are the ratios of the angular velocity of the 
field line to that of the BH horizon in the BZ and MC pro- 
cesses, respectively. 

Now we are going to discuss the relation between the 
power-law index S and the ratio ^out based on the transfer 
of angular momentum from region III to region II. 

It is assumed that the accretion rate (Atd) at rout is 

V / out 

related to (Ma) by 

('^'-)™.-(^^'')mc i'"^)siMC - Tbp-STmc, (40) 

where STmc = T^f(^ — Tj^^c ^^e angular momentum trans- 
ferred from region III to region II, and 5 is a fraction param- 
eter taken as 5 = 0.5 in calculations. The terms Tl"c and 
Ta{c ^-re the transfer rates of angular momentum at rsiMC 
and Tout, respectively. Substituting equation (15) into equa- 
tion (40), we have 

r,2 2 { S t 2 \ 12 



= Tbp - STmc- (41) 
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(a) (b) (c) 

Figure 4. The curves of Voutflow versus S for different values of a„ and n, (a) a* = 0.40, n = 7.6 (b) a* = 0.60, n = 6.6 and (c) 
a* = 0.80, n = 5.6. 




Figure 5. The curves of r]p versus S for different values of a* and n, (a) a* = 0.40, n = 7.6 (b) a* = 0.60, n = 6.6 and (c) a* = 0.80, 
n = 5.6. 
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Table 2. The values of the derived parameters with ^out/^siMC = 10^ 



a* 


n 


he 


^SIBZ 


S.SIMC 


S 


0.40 


7.6 


52.41 


12.67 


2.20 


0.0353 


0.50 


7.1 


41.55 


10.64 


2.21 


0.0355 


0.60 


6.6 


33.73 


9.31 


2.22 


0.0358 


0.70 


6.1 


27.40 


8.39 


2.23 


0.0362 


0.80 


5.6 


21.55 


7.74 


2.22 


0.0366 






(a) 



(b) 



(c) 



Figure 6. The inequality Pbp > Pbz with rjp > I holds in the shaded region for (a) S = 0.2, (b) S = 0.4 and (c) S = 0.6. 



Equation (|41|l can be rewritten as 



-.2 S 



" 'iSIMC 



Tbp - STmc- (42) 



Equations (33) — (42) provide a closed set for calculating 
the powers and torques in the BZ, MC and BP processes, 
in which seven parameters are a*, n, he, am S, S and ^out 
are involved. Besides the parameters Um = 0.1 and 5 = 0.5 
the rest five parameters are divided into two types. 

(1) The parameters a*, n and he are required to be 
greater than some critical values for the coexistence of SIBZ 
and SIMC as argued in W06, by which the boundary radii 
^siMC and ^sibz can be determined. 

(2) The parameters 5* and Cout are involved in the BP 
process, and only one of them is independent as shown in 
Figure 3 below. 

Thus we have four independent parameters (a*, n, he 
and (.out) as the input parameters, and ^simc, ^sibz and S 
as the derived parameters. Based on the above analysis we 
have the values of the derived parameters corresponding to 
the input ones as listed in Table 2. 

Based on equations (34), (39) and (42) we obtain the 
curves of the parameter S varying with Cout for the given 
values of o* and n as shown in Figure 3, and we find that S 
decreases monotonically with the increasing Cout- 

Incorporating equations (16) and (18), we have the ratio 
of Mout flow to Md expressed by 

1, (43) 



Tjoutflom = Mout flow/ (yMd) „ ~ 'io 



and the curves of rjoutfiow versus S for the given values of 
a, and n are shown in Figure 4. 



Inspecting Figure 4, we find that rioutfiow decreases 
monotonically with the increasing S, and this result implies 
that the greater ratio Cout corresponds to the less S and the 
greater outflow rate. 



4 JET POWER FROM BH ACCRETION DISCS 

As argued above the jets are driven respectively in the 
Poynting flux and hydromagnetic regimes by the BZ and 
BP processes, and the jet power can be estimated as the 
sum of the BZ and BP powers. 



Piet = Pbz + Pbp- 



(44) 



Thus we have the powers and torques in the BZ-MC-BP 
model as shown in Table 3. 

The results listed in Table 3 correspond to Cout ~ 10^, 
and we find that the BP power is greater than the BZ power, 
while the BP torque is about two orders greater than the BZ 
torque. Based on equations (33) — (39) we have the curves 
of rip = Pbp /Pbz versus S for the given values of a* and n 
as shown in Figure 5. 

Inspecting Figure 5, we find that the ratio rjp > 1 holds 
for a variety of values of the parameters, a,, n and S. Com- 
bining Figure 5 with Figure 3, we find that the BP power 
can be more stronger than the BZ power for less S and thus 
for greater Bout- 
in order to discuss the relative importance of the BP 
power relative to the BZ power in a visual way we plot the 
contour of = 1 in a, — n parameter space for the given 
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Table 3. The powers and torques in the BZ, MC and BP model with ^oiit 

= 102 



a* 


n 


he 


Tbz 


Tbp 


Pbz (xl0-») 


Pbp (xlO-^^) 


Pj.t (xlO-3) 


0.40 


7.6 


52.41 


0.077 


12.57 


4.00 


19.73 


23.73 


0.50 


7.1 


41.55 


0.089 


11.42 


5.94 


20.36 


26.30 


0.60 


6.6 


33.73 


0.098 


10.14 


8.18 


20.81 


28.99 


0.70 


6.1 


27.40 


0.107 


8.68 


10.89 


21.18 


32.07 


0.80 


5.6 


21.55 


0.118 


6.97 


14.69 


21.71 


36.40 



values of S as shown in Figure 6, in which the shaded regions 
indicate rjp > 1 for the BP power greater than the BZ power. 

Incorporating Figures 5 and 6, we find that the BP 
power could be greater than the BZ power for a wide range 
of the parametersa, , n and S. These results are consistent 
with the previous works in estimating the BZ power relative 
to the jet power from the inner disc (Ghosh & Abramowicz 
1997; Livio, Ogilvie & Pringle 1999). 

Cao & Rawlings (2004, hereafter CR04) estimated the 
jet powers of a sample of SCR PR I radio galaxies, and 
they argued that the BZ mechanism provides insufficient 
power to explain the high radio luminosities of at least a 
third, and perhaps all, of the sample, if the accretion discs 
in these sources are assumed to be advection dominated ac- 
cretion flows (ADAFs), or adiabatic inflow-outflow solution 
(ADIOS) flows. However, the BP power was not considered 
in CR04. From Table 3 we flnd that Pj^t = (2.5 ~ 6)Psz, 
and the jet powers can be amplifled significantly by using 
equation (44). 



5 DISCUSSION 

In this paper we incorporate three mechanisms into the BZ- 
MC-BP model, in which the jet can be powered by both 
BZ and BP mechanisms, which belong to the Poynting flux 
and the hydromagnetic regimes, respectively. In this model 
the ratio of the BP power to the BZ power varies with the 
concerned parameters, a*, n and S, as shown in Table 3 and 
Figure 5. 

As pointed out by Meier (2003) , there are observational 
reasons for believing that the same source may produce jets 
of rather different Lorentz factors, either simultaneously or 
in different accretion states. The results obtained in this 
model could interpret the jets of different Lorentz factors, 
and are consistent with a spine/sheath jet structure. The 
Poynting flux powered by the BZ process corresponds to 
the spine with higher Lorentz factor near the axis and hy- 
dromagnetic outflow powered by the BP process corresponds 
to the sheath with lower Lorentz factor away from the axis. 

Another feature of this model is the role of the MC 
process. It is assumed that disc accretion is depressed due 
to the transfer of the angular momentum from a rotating 
BH to the inner disc. It is the MC effects that give rise to 
the variation of the accretion rate as well as the outflow in 
region II. It has been argued that a very steep emissivity 
index required by broad Fe Ka lines can be interpreted by 
invoking the MC process, which is consistent with the XMM- 



Newton observation of the nearby bright Seyfert 1 galaxy 
MCG-6-30-15 (Wilms 2001; Li 2002c; W03). 

Recently, Sambruna et al. (2006) discussed the 
jet/accretion connection based on Chandra and XMM- 
Newton observations of three powerful radio-loud quasars, 
1136-135, 1150+497 (Chandra), and 0723-^679 {XMM- 
Newton), and they concluded that both jet and disc emis- 
sion contribute to the X-ray emission from the three quasar 
cores. It was pointed out that the beamed emission con- 
tributes roughly 50% to the total flux in 2-10 keV, while the 
disc emission dominates below 2 keV. We expect that the 
jet /accretion connection in these radio-loud quasars could 
be interpreted based on the BZ-MC-BP model, in which the 
jet power is equal to the sum of the BZ and BP powers, and 
the MC power plus the disc accretion power can be used 
to fit the disc luminosity. In addition, the different Lorentz 
factor in the jet can be fitted by adjusting the ratio of the 
BP power to the BZ power, and the broad Fe Ka lines can 
be simulated by invoking the MC process as argued in W03. 

Very recently, Wang, Ye & Huang, (2007, hereafter 
W07) fitted the twin peak high frequency quasi-periodic os- 
cillations (QPOs) associated with the jets observed in sev- 
eral sources, in which BH binary GRO J1655-40 is included. 
The fitting in W07 is focused on the association of the 3:2 
QPO pairs of the jets based on the coexistence of the BZ 
and MC mechanisms, in which the jet is powered by the BZ 
process and the MC process is invoked to fit the 3:2 QPO 
pairs. It is noted that an X-ray-absorbing wind discovered 
in an observation of GRO J 1655-40 must be powered by a 
magnetic process. Detailed spectral analysis and modeling 
of the wind shows that it can only be powered by pressure 
generated by magnetic viscosity internal to the disk or mag- 
netocentrifugal forces (Miller et al. 2006). Thus we expect 
that the fitting given in W07 should be improved based on 
the BZ-MC-BP model. 

Recently, a lot of works have been done on general rela- 
tivistic simulation on BH accretion and outflow (De Villiers, 
Hawley & Krohk 2003; Hirose, Krolik, De ViUiers & Haw- 
ley 2004; De Villiers, Hawley, Krolik & Hirose 2005; Kro- 
lik, Hawley, & Hirose 2005). McKinney & Narayan (2007) 
have found that a highly relativistic, Poynting-dominated 
funnel jet in the polar regions of a Kerr BH is associated 
with a strikingly simple angular-integrated toroidal current 
distribution dl^p/dr oc r"^''*, and the polar field is con- 
fined/coUimated by the corona. It is interesting to note that 
the BZ-MC-BP model is consistent with the simulations in 
some main features: Q Poynting-dominated jets driven by 
the BZ process, and ([2J a large-scale poloidal magnetic field 
varying with the disc radius as oc r^" corresponds to a 
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simple toroidal current with a power-law distribution in the 
disc. Wc shall improve this model based on the constraints 
of the numerical simulations as well as the observations in 
our future work. 
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